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ABSTRACT 

Compton scattering plays an important role in various astrophysical objects such as accreting black holes 
and neutron stars, pulsars, and relativistic jets, clusters of galaxies as well as the early Universe. In most of 
the calculations it is assumed that the electrons have isotropic angular distribution in some frame. However, 
there are situations where the anisotropy may be significant due to the bulk motions, or anisotropic cooling 
by synchrotron radiation, or anisotropic source of seed soft photons. We develop here an analytical theory 
of Compton scattering by anisotropic distribution of electrons that can simplify significantly the calculations. 
Assuming that the electron angular distribution can be represented by a second order polynomial over cosine 
of some angle (dipole and quadrupole anisotropy), we integrate the exact Klein-Nishina cross-section over 
the angles. Exact analytical and approximate formulae valid for any photon and electron energies are derived 
for the redistribution functions describing Compton scattering of photons with arbitrary angular distribution 
by anisotropic electrons. The analytical expressions for the corresponding photon scattering cross-section on 
such electrons as well as the mean energy of scattered photons, its dispersion and radiation pressure force are 
also derived. We applied the developed formalism to the accurate calculations of the thermal and kinematic 
Sunyaev-Zeldovich effects for arbitrary electron distributions. 

Subject headings: accretion, accretion disks - cosmic background radiation - galaxies: jets - methods: analyt- 
ical - radiation mechanisms: nonthermal - scattering 



1. INTRODUCTION 

Compton scattering is one of the most important radia- 
tive process that shapes the spectra of various sources: black 
holes and neutron stars in X-ray binaries, pulsars and pulsar 
wind nebulae, jets from active galactic nuclei, and the early 
Universe. Compton scattering kernel takes a simple form if 
electrons are ultra-relativisti c with the Lorentz factor y » 1 
(iBlumenthal & Gouldlll97C)h . In a general case, when no re- 
stricti ons are made on the energies of photon and electrons, 
I Jones! (1968) derived the kernel for isotropic electrons and 
photons. The formulae there co ntain a fe w misprints, but even 
by correcting those (see e.g. IPe'er & Waxman 2005) they 
cannot be used for calculations because of a number of nu- 
merical cancellations (see e.g. iBelm ont 2009). An alternative 
derivation to that kernel was given by Brinkmann (119841) and 
iNagirner & PoutanenI ( 11994 NP94 hereafter), who showed 
how to extend numerical scheme to cover all photon and elec- 
tron energies of interest in astrophysical sources. 

In real astrophysical environments, the radiation field does 
not need to be isotropic and a more general redistribution 
function is required to describe angle-dependent Compton 
scattering. NP94 extended previous results to the situation 
when the photon distribution can be represented as a linear 
function of some polar angle cosine (Eddington approxima- 
tion), derivi ng an analytical formula for the first moment of 
the kernel. lAharonian & Atoyaiil (Il981h were first to derive 
a redistributio n function for arbit r ary photon angular d epen- 
dence (see also ' Prasad et al.lll986l) . iKershaw et"al] ([1986) and 
iKershaw ( 1987) have developed numerical methods to com- 
pute the kernel efficiently and with a high accuracy. All these 
w orks neglect the effect of pho ton polarization. 

INagirner & PoutanenI (119931) have derived a general Comp- 
ton scattering redistribution matrix for Stokes parameters as- 
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suming an isotropic electron distribution. A general rel- 
ativistic kinetic equation incorporating the effects of in- 
duced scattering and polarization of photons as well as 
electro n polarization and d egeneracy has been derived by 
Nagir ner & PoutanenI (120011) . 

In this paper, we propose a method to extend previous re- 
sults to the case where the electron distribution is no longer 
isotropic, but can have weak anisotropics which can be repre- 
sented as a second order polynomial of the cosine of some 
polar angle. The proposed formalism can find its applica- 
tion in a number of astrophysical problems. The distortion 
of the cosmic microwave background (CMB) caused by the 
hot electron gas in clusters of galaxies (i.e. kinematic and 
thermal Sunyaev-Zeldovich effect) is an obvious application. 
The electron distribution, isotropic in the cluster frame, can 
be Lorentz transformed to the CMB frame, resulting in a 
dipole term linear in cluster velocity and a small quadrupole 
correction. Compton scattering then can be directly com- 
puted in the CMB frame. Another possible application con- 
cerns th e models of o utflowing accretion disk-coronae or jets 
(iBeloborodov 1999:,Malzac et al. 2001) . If the outflow veloc- 
ities are non-re lativistic, the radiation transport can be consid- 
ered directly in the ac c retion disk frame following recipe by 
iPoutanen & SvenssonI (119961) . with the Lorentz transformed 
electron distribution. 

Weak anisotropics in the electron distribution can 
arise in high-energy sources with ordered magnetic field, 
because of the pitch angle-dependence of the syn- 
chrotron cooling rate and/or anisotropic injection of high- 
energy electrons (e.g. Biornsson 1985t Roland et al.' 119851 : 
Crusius-Waetzel & Lesch 1998; Schopper et al. 1998). An 
anisotropic electron distribution is also a very natural outcome 
of the photon breeding operation in relativistic jets wit h the 
toroidal magnetic field (IStern & PoutanenI 120061 120081) . be- 
cause the electron-positron pairs born inside the jet by the ex- 
ternal high-energy photons move perpendicularly to the field. 
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Our method is also extendabl e to the polarized radia t ion us - 
ing the techniques developed in lNagirner & Poutanenl(ll993h . 
It is also in principle possible to calculate the scattering redis- 
tribution function in the case when the electron distribution is 
expressible as an arbitrary order expansion over the polar an- 
gle cosine. Unfortunately, in the latter case the analytical ex- 
pressions become extremely cumbersome and the advantage 
over direct numerical integration becomes small. 

Although here we consider only photon scattering it is also 
possible to apply the same method for electrons interacting 
with the photons in the case when photon angular distribution 
is expressible as an expansion of powers of the polar angle 
cosine. This can be of interest only in the deep Klein-Nishina 
regime where the electron can lose or gain a large fraction of 
its initial energy in one scattering and continuous energy loss 
approach is not applicable. 

Our paper is organized as follows. In Section |2] we intro- 
duce the relativistic kinetic equation for Compton scattering 
and define the redistribution function and total cross-section. 
The expressions for the total cross-section, the mean energy 
and dispersion of scattered photons, and the radiation pressure 
force are given in Section |3] The exact analytical formulae 
for the redistribution function for mono-energetic anisotropic 
electrons as well as approximate formulae valid in Thomson 
regime are derived in SectionH] We present the numerical ex- 
amples of redistribution functions in Section|5] where we also 
develop the relativistic theory of Sunyaev-Zeldovich effect. 
We summarize our findings in Section|6] 

2. RELATIVISTIC KINETIC EQUATION 

Let us define the dimensionless photon four-momentum as 
X = {x,x} - x{l,co}, where cj is the unit vector in the pho- 
ton propagation direction and x = hvjm^c^. The photon 
distribution will be described by the occupation number n. 
The dimensionless electron four-momentum is p = {y, p] - 
{■y,pfl} = y{l,y6Q), where Q is the unit vector along the elec- 
tron momentum, y and p - -y/y^ - 1 are the electron Lorentz 
factor and its momentum in units of m^c, and /3 - p/j is the 
electron velocity in units of speed of light. The momentum 
distribution of electrons is described by the relativistically 
invariant distribution function fe{p) (see iBelvaev & Budken 
[1956; NP94). 

The interaction between photons and electrons via Comp- 
ton scattering (in linear approximation, i.e. ignoring induced 
scattering and electron degeneracy) can be described by the 
explicitly covariant relativistic kinetic equati on for photons 
dPomraningll 1 9731: iNagirner & Poutanenll 1 9931: NP94): 
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where V - {d/cdt,V} is the four-gradient, re is the clas- 
sical electron radius, F is the Klein-Nishina reaction rate 
(iBerest etskii et al. 1982) 
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are the four-products of corresponding momenta. The second 
equahties in both equations © arise from the four-momentum 



conservation law. The invariant scalar product of the photon 
four-momenta can be written in the laboratory frame as well 
as in the frame related to a specific electron 



q = X-X^ = XXi(l -yu) = ^^i(l -jUo) = ^-^1, 



(4) 



where p - oj-Mi is the cosine of the photon scattering angle in 
some frame and fiQ is the corresponding cosine in the electron 
rest frame. 

In any frame, the kinetic equation can be also put 
in the usual form of t he radiative transfer equation 
(iNagirner & Poutanenlll993l) : 

1 <5 \ 
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c ot I 
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where A^e is the electron density in that frame. Here we have 
defined the photon redistribution function 
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and the total scattering cross-section (in units of Thomson 
cross-section crx) 

'd^pd^px d^xi rr ^ X 
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2.1. Electron distribution and scattering geometry 

Let us now consider a specific frame E. Our basic assump- 
tion is that the anisotropy of the electron distribution in this 
frame can be expressed as a second order polynomial expan- 
sion in the cosine of the polar angle in some coordinate system 



1 ^ 

" k=Q 
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where A^e is the electron density in frame E, rj^ - £1 ■ I3 
is the cosine of the polar angle of the electron momentum, 
Pki^le) are the Legendre polynomials, and we now switched to 
the dimensionless distribution function /e(y, //e) normalized to 
unity: 



/e(%?7e) P dp= 1. 



(9) 



The moments /o, /i and /2 determine the energy spectrum 
of electrons and the relative magnitudes of the isotropic and 
anisotropic components. The distribution function for mono- 
energetic electrons of energy yo can be described by 
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where the ratios /1//0 and /2//0 are constants. 

The directions of photons in this coordinate system (see 
Fig.[T]) are given by 

oj- -yjl - 77- COS0 /i + 1 - 77- sin /2 -H 77 /3 , (11) 

'^i = -y/l -?7I COS 01 li + -^1 - T]l sin0i I2 + h, (12) 
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Fig. 1. — The scattering particles' momentum vectors and the vector /3. Note 
that the shown angular variables are the cosines of the respective angles. 

SO that the cosine of the scattering angle is 

^Ji ^ <o ■ (jJx ^ rjr]i + - if ^\ -rj]cos(4>-4>i). (13) 

3. TOTAL CROSS SECTION AND MEAN POWERS OF ENERGffiS 

3.1. Total cross-section 

Let us simplify the expression for the total cross-section. 
We follow here the approach described in NP94. We rewrite 
the cross-section as 



so(x) 

where 

Using the identity 



5o(^)^/e(r,'7e)— , (14) 
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6{y\ + xi - y - x) ^ ji S[x^ ■ {p_ + x) - X- p^ 
and taking the integral over pi in equation ( fTSl l, we get 
3 1 fd^xi / V 



3 



Fd^i 



(17) 



f/(l+2f) 



where we used invariant x\dx\d^(jj\ - ^\d^\ djiQ d(pQ 
and the fact that F does not depend on azimuthal an- 
gle (1)0. Substituting F from Equ ation (|2]i we get 
(iBerestetskii. Lifshitz. & Pitaevskiilll982t NP94) 
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(18) 

Putting ^ ^ X, we, of course, get the total Klein-Nishina 
cross-section for a photon of energy x on electrons at rest. 

To obtain the total scattering cross-section on an anisotropic 
electron distribution, we have to calculate the angular inte- 
grals over incoming electron directions in Equation (fT4] i. We 



introduce the cosines between electron momentum and pho- 
tons: 

^=Q-6J, (19) 

so that 



^^x{y-pC), ^1 = xi(y-p^i). 



(20) 



We choose the spherical coordinate system and measure the 
polar angle from the direction of the initial photon cj, we get 



In 



MX)^MX, '?) = ^ J ^ J H SoiO ^JdO /e(r, /?e). (21) 

-1 

Azimuth C) is now defined as the difference between the az- 
imuths of the electron momentum direction and the vector /3 
in a frame with z-axis along a>. The angular variable rjs in the 
expansion ([8]l is expressed in this frame as 



77e = 7?^ + ^l-lf ./W^COS O. 



(22) 



The physical meaning of Equation ( 1211 1 can be also under- 
stood if we consider a monoenergetic beam of electrons along 
/3 axis: /e(y, 7/e) = ^(T/e - l)^(r - ro)/(27r/?r). Then 



soix^rf) = (1 - prf) sq(x'). 



(23) 



where x' - xy{l -prj) is the photon energy in the electron rest 
frame (we omitted subscript in y and fi). The factor I - firj 
in Equation ( |23] | accounts for the reduced number of collision 
per unit length. 

When calculating the azimuthal integral in equation (1211 1 
we just have to integrate Pkirje) with ?/e given by Equation 
( |22] |. The properties of the Legendre polynomials give us the 
average 

KUh) = Pk(r]) PdO, (24) 
so that the averaged distribution function becomes 

1(7,^1,0^^ fd<i>My,7j,)=yMy)Ptiv)Pi<(0- (25) 
The cross-section now can be written as 
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Xoo 
py dy fk Aok(x,y), (26) 
K=y, 



where 



1 r' 

^ok(x, y) = - — 

2yx J_i 



PkiO^^oiOdi. (27) 



Changing the integration variable and expressing Pk through 
^ using Equation ( l20t . we get 



Aoi:(x,r) = ^ bnkXOn, 



(28) 



«=0 



where 



Xon(x,y) = —^f' r'^o(^)d^. (29) 



and 
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Fig. 2. — (a) Total cross-section XKs = Aqo for isotropic mono-energetic 
electrons of various momenta p = 0.1, 1, 10, 10", lO' (from top to bottom, 
dot-dashed, solid, dotted, dashed, and triple-dot-dashed curves) as a function 
of photon energy x. (b) Relative correction to the cross-section arising from 
the dipole term in the electron distribution (s) for the same p as in panel (a). 
At small X the curves approach the asymptotic value in Thomson limit -j8/3. 
(c) Relative correction to the cross-section arising from the quadrupole term 
A02/A00. 



The zeroth function Aqo = Xoo coincides with the function 
^oix, 7) from NP94. When electron distribution is isotropic 
(i.e. f\ - fi - 0), expression (l26b for the total cross-section 
is reduced to equation (3.4. 1) from NP94 and the dependence 
on T] obviously disappears. Functions xon of two arguments 
can be presented through the functions of one argument; 



Aron(jc,r) 



1 U 



2+n 



where 



2jp 2 + n 
i+1 



iA«+i,o(Jc«) 



t'so(t)dt. 



(31) 



(32) 



The explicit expressions for the functions i/',()(^) (/ = 1,2,3) 
can be found in Appendix iB] (see also NP94). Thus the to- 
tal cross-section is given by a single integral over the electron 
energy ( l26b with all functions under the integral given by an- 
alytical expressions. Numerical calculations of functions Aqa 
can be separated into three regimes: (1) in Thomson regime, 
jcy <K 1, the series expansion (see Appendix iDb can be used; 
(2) for <K 1, but X not sufficiently small for regime (1), we 
numerically take the integral in Equation dZTl i using 5-point 
Gaussian quadrature to reach accuracy better than 1%; (3) in 
other cases, we use the sum in Equation ( |28] l and analytical 
expressions dSTT i for;t'o«. 

For mono-energetic electron distribution of Lorentz factor 
7 given by Equation ( fTOb . we can introduce the cross-section 
analogously to Equation (|26] |: 



fk 

~)(x,y,T]) YPk(T])AQk(x,y). 



(33) 



For isotropic mono-energetic electrons, the total cross-section 
is shown in Figure |2^. The relative corrections arising due to 
the dipole and quadrupole term in the electron distribution are 
shown in Figures |2}3 and |2];, respectively. These have to be 
multiplied by the angle- and, possibly, the energy-dependent 
factor Pkiififkl fo to obtain the final correction. In the Thom- 
son limit, at small xy <^ 1, the cross-section takes the form 
(see AppendixiDli 



SQ{x,y,rf} ~l-~r]/3, 
3 /{) 



(34) 



where the correction to unity term can be easily obtained by 
averaging the transport cross-section over electron directions 
(i.e. integrating (1 - POrj^ over the angles). This corresponds 
to the flattening in Figure |2]3 at Aoi/Aqo = -j0/3. The cor- 
rection from the quadrupole term in this regime as well as for 
non-relativistic electrons becomes negligible: 



Ao2/A()o ~ -y^p (-^r)- 



(35) 



3.2. Mean powers of scattered photon energy 

In some situations, the full relativistic kinetic equations can 
be substituted by the approximate one obtained in Fokker- 
Planck approximation. This requires knowledge of various 
moments of the redistribution function, such as total cross- 
section, the mea n energy and dispersion o f the scattered pho- 
tons (see NP94, 1 Vurm & Poutanenll2009l) . It is often time- 
consuming to compute numerically the integrals of the redis- 
tribution function and instead direct calculations of the mo- 
ments are preferable. Below we obtain analytical expressions 
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for the mean energy and dispersion of the energy of scattered 
photons in frame £ as a function of the initial photon energy 
X and the direction of its propagation relative to a symmetry 
axis of the electron distribution h. 

Following NP94, we define the mean of powers of energy 
of scattered photons: 



}so(^^M7,il.)—, (36) 

r 



where now 



<^{)*o© = i x{ 5\p^ +x,-p-x) (37) 

167T^J yi xi i -1 1 - 



3 1 

167?^ 



■■'^^dxicfui F S{xi[y + x - (Oi -(pQ+xw)] - ^) 



3.2.1. Averaging over photon directions 

Quantities (ITtT i are not invariants (except for j - 0), and 
we have to compute the scattered photon energy is a certain 
frame, which we choose to be frame E. Because of the addi- 
tional term Xj under the integral, a simple change of variables 
to the electron rest frame as in Equation ST% is not possible. 
Instead, we use the 5-function to take the integral over xi: 



-^1 fx 



■f 2 d^w, F . 



(38) 



Now we change the variables to those in the electron rest 
frame (with subscript 0). We choose the coordinate system 
with the polar axis along the direction of the incoming photon 
a>{). In this frame, the cosine of the angle between the elec- 
tron momentum and the incoming photon is ^o- The cosine 
of the angle between the outgoing photon momentum and the 

electron is then = ^o/Uo + -^1 - ^1 - f4,cos 0o- 

We use invariants Xjd^iwi - ^^d//()d0o and the energy con- 
servation law in the electron rest frame = ^/(l H-f[l -fio]), 
to get (see NP94) 



Xid^wi - d^id^o- 



Finally, we have 



1671 f Jf/(I+2f) Jo 



(39) 



(40) 



Because is the energy of scattered photon in the electron 
rest frame, the Doppler effect gives us 



(41) 

where we now can substitute 



, 1 1 

UO = 1 H 

^ ^1 



P^Q " I ~ ^' 



(42) 



which are consequences of the conservation law and of the 
Lorentz transformation x = ^{y+ p^o), respectively. The terms 
containing a linear combination of square roots and cos 00 
will disappear after averaging over (po. 
We now introduce moments of the invariant cross-section 



f/(l+2f) 



(43) 



For j - 0, we get of course the total cross-section sq given by 
Equation ( fTSl ). NP94 derived the corresponding expressions 
for = 1,2: 



(44) 
(45) 



where = In (1 -n 2 ^), and /Jf = 1/(1 + 2 ^). 
For the mean energy of the scattered photon we have then 

{xi}so(^)^^[ySi(0 + xS2m, (46) 
and for the mean square of energy 

<x2) ^o(^) = y'f S4(0 - yx^S^iO + x^ Se(0 - f S^^), 

where 



(47) 



5 1(^) = [sq(^) - Slim S2(^) = [*i(^) - Si(£)] 

Si(^) = [sd^) - S2(^)] /f, Sd^) = [5 1(^) - 53(^)] /^, 

Ss(^ = 3 S4(^) -4S3(a SA^) = S,(^) - S 4(^)12, 
SdS)^s2{^)-2>S,{^). (48) 

All functions S are elementary. In addition, they are de- 
fined in such a way so that not to become zero at ^ = 0. The 
series expansion of functions s and S for small arguments are 
presented in AppendixlAl 

3.2.2. Averaging over electron directions 

We need to integrate in Equation ( |36] | over anisotropic elec- 
tron distribution. We follow the derivation of the total cross- 
section that lead from Equation (ISTT i to Equation (l26T l. Repre- 
senting integral over electron momentum d-^;? = p^d/^d^dd) 
we get: 



^ f\OS 

x( i^(x, 77) = 47r X-' ^ Pkirf) I py dy ft Ajk(x, y) 
where 

If' * 

A>i(x, r) = 5^ J ^ PkiO ^ iA> ^o(^) = g 

and 



(49) 



nk Xj n 

(50) 



Xjn{x,y) = 



2ypx 



2+j+n 



L 



■x{y+p) 



{x{)s^U) df- (51) 



Functions ;t'jO - A^o coincide with functions ^^^(x, y) in- 
troduced by NP94, while functions X(in are given by Equa- 
tion OH). The explicit expressions for the function Xjn for 
j - 1,2 (which are analogous to functions '^i and *i'2 from 
NP94) can be obtained using expression for mean powers of 
energies (|46] | or (l47b : 



1 



,3+n 



X\n{x, y) = [y'^2+n,\ixu) + X^2^n,2ixu)\ 

lyp i + n 



(52) 



u=y-p 



X2n(x,y) = 



1 

2yp 
2 + n 



,,4+n 



,3+n 



— 'V3+n,4(xu) - y- *I'2+«,5(JCm) 

■ n i + n 



"Vi+nfiixu) ■ 



.A+n 



A + n 



^3+nj(xu) 



u=y+p 



(53) 
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Fig. 3. — (a) Mean energy of scattered photons in units of the incident pho- 
ton energy l(i"'°(x,y)/x = A|o/Aoo as a function of x for isotropic mono- 
energetic electrons of various momenta p = 0. 1, 1, 10, 10^, lO' (from bottom 
to top, dot-dashed, soUd, dotted, dashed, and triple-dot-dashed curves). The 
asymptotic value at small x in Thomson approximation is 1 4- 1 . (b) A cor- 
rection to the mean energy (in units of Tf "'°) arising from the linear term in 
the electron distribution (8) with fi/fo = 1 . SoHd, dotted and dashed curves 
correspond for p = 1, 10, 100, respectively. The curves from top to bottom 
con'espond to = -1, -0.5, 0, 0.5, 1. At small x, the curves approach the lim- 
iting value given by Equation (58). (c) Same as (b), but for the quadrupole 
term in the electron distribution (8) with /2//0 = 1. These are even functions 
of the curves from the bottom to the top correspond to f; = 0, 0.5 , 1. At 
small X, the curves approach the limiting value given by Equation (59). 



where 

These are related to functions 



Jo 



(t)dt. 



(54) 



(55) 



because functions S j are expressed through sj. The explicit 
expressions for both type of these functions as well as their 
series expansions for small arguments are given in Appen- 
dices |B] 

As in the case of functions Aok, for calculating Ajk, we con- 
sider three regimes: (I) xj <s: 1, when we use the series ex- 
pansion (see Appendix|D]i; (2) for <K 1 we numerically take 
the integral in Equation ( ISOl l using Gaussian quadrature; (3) 
in other cases, we use the sum in Equation ( fSOl l and analytical 
expressions for xjn- 

For mono-energetic electron distribution (fTOl i of Lorentz 
factor y, we can introduce the mean powers of photon energy 
analogously to Equation (|49] i: 



X\ So(x, 7, T]) 



x^Y^jPM^. 



■jk- 



(56) 



The mean energy of scattered photons for such electrons for 
a scattering act is given by the ratio of Eqs. ( |56] | and ( |33] ). 
It is shown in Figure [3^. In the low-energy (Thomson) limit 
the energy gain factor is given by a well known expression 
'xi"'°/x = I + 4-p^/3, which translated to |y- at large y. The 
relative corrections arising due to the dipole and quadrupole 
term in the electron distribution are shown in Figures [S]^ and 
[3]:, respectively. Using asymptotic expansions of Ajk in the 
Thomson limit (see AppendixlDl). we get the asymptotic value 



Xl(x, J, T]) _ T 
— iso/ \ 

xi (x, y) 



(l + i/?2)-|^y2A;,+ ^/?Vff2(/7) 



r 



(57) 

Thus in non-relativistic limit /? <K 1, the correction is negli- 
gible. In the relativistic limit 7 » 1, the relative corrections 
arising from the two terms are 



xi(x, y, rj) 

ir"°(x,r)" 1 

xlix, y, T]) 



1 - 



2fo' 



(x,y) 10/0 



(58) 
(59) 



3.3. Energy exchange and dispersion 

The difference of the photon energies before and after scat- 
tering x-xi is of course just the energy transfer to the electron 
gas. For the fixed angle between electrons and incident pho- 
tons ^ (and fixed electron energy y), the energy loss averaged 
over the directions of scattered photons is x - (xi). The prod- 
uct {x-{xi))NgCrj i()(^) is then the energy loss on a unit length. 
From Equation ( l46t we can easily get (see also NP94): 

(x-{xi))so(^)^xso(0-y^Si(0-x^S2(^) = (x+x^-y^)S 

(60) 

The corresponding energy loss (per unit length and in units 
NeO-j) averaged over the electron directions (and integrated 



Theory of Compton scattering by anisotropic electrons 



7 



o 




o 



0.2 



0.1 



^ 0.0 

w7 



-0.2 



: ' 1 ' 1 ' 1 
1 (b) 

: / \ L. 


1 : 

-■' '■. / ^ 
'. ' * 
■-. / ' 

\ .' ..V * - 
\ .-' ■■ / ■. ' ^ - 
.■■\/ V ^ ^ - 
\ \ / '.''■ ^ \ ~ 


' ^ ^'"^y 


// 

7 '■ 1 1 - 
I ■ .- \ ■ \ , , - 

/'■• ■■■■ / f I 

V = 


//' 'x .■'/'/ y \> 


~ #■ / A 

■'"'n ■■' / ■ \ 


I~ ' '\ ■■ / ' \ 

' .^ y / ■ \ y 





10-' 



lO"" 



10"^ 



10" 10^ 



0.1 



' 1 


'I'I 


' 1 






(c) : 




^ '- *\ 
^ '• ' \ 




^ ■•' ' 1 \ / 

/ i / Vy 
^ ' i / 




1 


\ , 1 


1 



10' 



10-* 



10^ 



10" 



10^ 



Fig. 4. — (a) Dispersion of the energy of scattered photons (in units of jr^) 
for isoti'opic mono-energetic electrons of p = 0.1, 1, 10, 10^, lO' (from bot- 
tom to top). The asymptotic value of Dj^ at small x in Thomson approx- 
imation is ^(237^ - 8)p^. (b) A correction to the dispersion (in terms of 
the isotropic quantity) arising from the linear term in the electron distribu- 
tion (s) with /i//o = 1. Solid, dotted and dashed curves correspond for 
p = 1, 10, 100, respectively. The curves from top to bottom correspond to 
= -1,-0.5,0,0.5,1. (c) A correction to the dispersion arising from the 
quadrupole term in the electron distribution (s) with /2 //o = 1 for the same p 
and as in panel (b). These are even functions of the curves from bottom 
to top coiTespond to f; = 0, 0.5, 1. 



over electron energies) becomes [see Eqs. ( |26] | and (|49]l1: 

(x - xD 'so(x, 77) = 4 TTx '^Pkiri) I py dy fk (Aq^- - Ah) . 

(61) 

The heating rate per unit volume is then 

E = A^eo-T ^ ^ I(x, w) ^1 - ^ j ^{x, 77), (62) 

where I{x, a>) is the specific intensity of radiation in a given 
direction. This expression can be positive (so called Comp- 
ton heating) when the photons typically have larger energies 
than the electron gas, or negative (Compton cooling) when 
one considers cooling of the relativistic electron gas by soft 
radiation. 

The dispersion of the scattered photon energy is given by 

the usual expression D(x) - Xj -Icf', which of course depends 
on the electron momentum distribution. For mono-energetic 
electrons we can define the dispersion as 



D(x,y,T]) ^ x1:(x,y,T]) - xi {x,y,r]). 



(63) 



where Xj(x, 7, //) are given by Equation (I56b . The dispersion 
for isotropic electrons is shown in Figure [3^. The low-energy 
(Thomson) Umit for x <s: 1/y is (see NP94) 



D(x,y)^x'- (23 r'- 8) 



(64) 



The relative corrections arising due to the dipole and 
quadrupole term in the electron distribution reach about 50 
per cent and are shown in Figures andH):, respectively. 

3.4. Radiation force 

Now we would like to derive analytic expression for the 
radiation force acting on the electron gas. NP94 have devel- 
oped a formalism appropriate for isotropic electron distribu- 
tion, when the averaged transferred momentum is along the 
momentum of the incoming photons, because of the symme- 
try. For the electron distribution described by Equation 
the momentum is transferred in the plane containing the initial 
photon momentum and the symmetry axis I3. If the incident 
photons are axially symmetric around I3, then obviously, the 
total momentum transferred to the electrons has be parallel to 
/3 by symmetry. We derive here more general formulae for 
the total momentum transferred by a beam of photons prop- 
agating along direction ai (such as a> ■ I3 - rj), as well as its 
projections to /3 and perpendicular direction. 

Let us introduce the vector basis: 

It, — rioj (i) X It, 

ei(cj,h) = , e2(ci),h) = , C3 = oj. (65) 

In a single scattering act the momentum transferred is 

Q - xoj - xi(Oi. (66) 

The components of the momentum transferred to the electron 
gas along and perpendicular to o) are: 



Qt, — X — xi oji ■ ej, — X — xi/j., 

m -w 



Q 



-xi oj\ ■ e\ 



-x\ 



(67) 
(68) 
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Analogously to Equation ( |36] ), we define the mean transferred 
momentum as 



Poutanen & Vurm 

and using identity = (jx - ^)/p, we finally obtain 



(e)^0(^)^/e(r,'7e)— , 

7 



(69) 



where 



1671^ J yi 



d^jc, 



_ 3 1 r 



xidxid cji FQ dixilj + x - (Oi-{p£l+x(o)] - ^} . 



3.4.1. Averaging over photon directions 

Let us introduce a vector basis defined by the photon and 
electron momenta: 



e\{(o, SI) - 



^3 = 6>, (71) 



where - o) ■ Q., therefore n = -y/l - ei((y, Jl) + ^03. 
Fixing the angle arccos ^ between the electrons of momentum 
p il and the incident photon momentum and averaging over 
directions of scattered photons, we can get the mean momen- 
tum transmitted in the direction a>: 

{Q3} s^M) = <x - xi}so(0 + {xi (l-fi)} io(^). (72) 

The first term is given by Equation ( l60b , the second term is 

3 



lOTTX^ J 



(1 -;u)Fd^wi 



(73) 



_ 3 

" 8^ id 



(f-^i)Fd^i = -[so(^)-sd^)] = ^5i(a 

f/(l+2f) ^ 

where we have used the invariant given by Equation (01) and 
changed the variables according to Equation ( [39] ). Thus, for 
the fixed electron and photon energies and the angle between 
their momenta, the mean momentum transmitted to the elec- 
tron gas in the direction of the initial photon propagation o), 
in accordance with (gHJ, (|46ll and (|73]l, is (NP94) 



x + x^-r^+- 5i(a (74) 

X ' 



In contrast to NP94, we are now interested to know the to- 
tal momentum transfer Obviously, by symmetry, it has to lie 
in the (£1, o)) plane. Averaging expression ( 1681 ) for Qi over 
angles is not easy, but we can compute the momentum trans- 
ferred along the electron momentum: Qn = xl^ - x\t^\. Using 
identity x\l^\ - (yx\ - ^x)jp and substituting Eqs. ( 1411 1 and 
(l42l) . similarly to Equation ( |40] |. we get 



3 1 r 

<xi^i)^o(f)= — - 
3 



f/(l+2f) 



X 

Jo 



1 

Fdfi- 



yx 
J 



yx 



I— + ;^---iKi+r|r-7 



1 



= - [y^^SiiO + yx^S2(0-^si(0] , 



(75) 



<en).^o(^) = {x^-xi(i)so(^) ^-(yx + yx^- y^^ - f) 5 i(f). 

(76) 

The momentum along e](oj, £2) is then simply 

/r, \ tt:\ ^2n> - i {Qi) 

{Q\)si.M) = , — soiO 



l^L-2yL + i.\^-pJr^2^Sm- (77) 



p^fT^C 

Thus the total transferred momentum Q averaged over direc- 
tions of scattered photons can be decomposed into two com- 
ponents along basis vectors: 



3.4.2. Averaging over electron directions 



(78) 



As in previous sections, we choose to measure azimuth of 
the electron momentum O in the frame defined by Eqs. ( l65T l 
from the projection of vector It, onto the plane perpendicular 
to 10. Therefore, 

<e> = <ei>[cosOei(w,/3)-HsinOe2(w,/3)]+<e3)e3. (79) 

The total momentum transfer avera ged over the electron dis- 
tribution is obtained from Equation ^69]): 

1 r'^TT 

e^(x,77)=- pdy d<^{Q}so(OMy,T],). 
X Ji J-i Jo 

(80) 

Obviously, the component along e2i(o, I3) becomes zero, as /e 
is an even function of The term along e 1 (o), Ix) involves 
integration of /e over azimuth with the weight cos O, and its 
averaged value is 



where P| are the associated Legendre functions: 



P\iu) = Vl -m2, pUu) = 3 m Vl -m2. 



(81) 



(82) 



It is worth mentioning that the isotropic component of the 
electron distribution /o does not contribute to the momentum 
transferred perpendicular to (o by symmetry. Substituting ex- 
pression (iTTl i to Equation (I8OI 1. we get the first component of 
vector 

QMx,r])^^7TxY^Pl(jj) I pydyfkAf,{x,y), (83) 



where 

^fk(x,y) 



1 



j_f' 

yx^ J-i 



Pl(0^iQi}^o(Od(. (84) 



k(k+l)2yx^ j_i 
Changing the integration variable to ^, and introducing a set 
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of functions 



1 n->:(y+p) I 



1 1 

2ypp 
5 + n 



•Xr-p) 

^2+n,\(xu) - 2y- *I'3+„j(xm) 

3 + n A + n 



u=y+p 



"=y-p 



« = 0, 1, 



we get 



(85) 



(86) 



Now let us evaluate the component of vector dSOl l along ct,. 
Because {Q^) does not depend on azimuth the azimuthal 
integration just gives the averaged electron distribution given 
by Equation dZSl l. Thus we get 



Qi soix, t])^4jtx ^ Pk(r]) j py dyfk (a, 



■Ok 



■ Au + 



(87) 



where 



1 I. 

1 r* 

Au(^,r) = :r^ M'i(0^<^i(l -M))^o(^)d^=y fe„i;ti;, 



and 



(88) 



1 nxir+p) 

^■"^"'^^^2^;^ <xi(i -^)).o(f)r' d^ 

^7P X Jx(y-p) 



1 

■ *J'3+« i('''^m) 

2yp 4 + n 



u=y+p 



u=y-p 



n = 0,1,2. (89) 



For isotropic electron distribution, the only function of in- 
terest is Ajg which coincides with function *F*(x, y) intro- 
duced by NP94. Now combining Eqs. dSST l and (IHtT i. we get 
the momentum transfer along the symmetry axis of the elec- 
tron distribution It, and perpendicular to it: 



Q±=-^Qi + V^^a- 



(90) 
(91) 



Expressions ( l87b . ( 1831 ) and ( |90l l give the momentum trans- 
ferred to the electron gas (in terms of one integral over the 
electron energy) along m, perpendicular to that direction as 
well as along vector Ij and perpendicular to it. 

Similarly to Equation ( |62| |. we can also get the two compo- 
nents of the momentum transfer rate per unit volume: 



PU3 = 



d^ojI(x,i])Qi^,Mx,T]). (92) 



For mono-energetic electron distribution of Lorentz factor 
y given by Equation (fTOl i. the momenta transferred along w 




Fig. 5. — (a) Average momentum transferred to the electron gas per scat- 
tering Q (in units of a) for isotropic mono-energetic electrons (with /i = 
/i = 0). Curves from bottom to top con'espond to the electron momenta 
p = 0.1, 1, 10, 10^, 10^^. The asymptotic value at small x in Thomson ap- 
proximation is 1 + 2p^/3 [NP94; see Equation )95H . (b) The momentum 
transferred along b> for anisotropic electrons with /i//o = 1 in units of 

the isotropic quantity Q . Solid, dotted and dashed curves correspond to 
p = 1, 10, 100, respectively. The curves from top to bottom correspond to 
= -1,-0.5, 0, 0.5, 1 (pink, green, black, red, and blue curves, respectively), 
(c) Same as (b), but for the electron distribution {8) with the quadrupole term 
with /2//0 = 1 for the same p and as in panel (b). These are even functions 
of the curves from the bottom to the top correspond to = 0, 0.5, 1 . T he fla t 
parts of the curves correspond to the Thomson limit given by Equation )95t . 
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10° 10^ 



Fig. 6. — (a) Momentum transferred in the direction perpendicular to o) (in 
units of tlie momentum along d for isotropic distribution) arising from the 
linear term in the electron distribution {8) with /i//o = 1- For rj = ±\, the 
momentum is zero by symmetry. From the top to the bottom curves corre- 
spond to = -0.5 (green), 0.5 (red), (black). Solid, dotted and dashed 
curves correspond to p = 1,10, 100, respectively, (b) Momentum transferred 
in the direction perpendicular to o) arising from the quadrupole term in the 
electron distribution (8) with /2//0 = 1 for the same p as in panel (a). The 
curves from bottom to top correspond to f; = ±0.25, 0.5, 0.75 (black, red, 
blue curves). The momentum is zero for = -1,0, 1 because of the sym- 
metry. Th e fla t parts of the curves coiTespond to the Thomson limit given by 
Equation (96). 



and perpendicular to it are 

2 



i=0 

2 



a ^(x, j,j^)^xj]f Pkiii) (a,u - Ai, + a;,) , (93) 

1— n ■'0 



k=\ 



(94) 



where we kept the notations for the functions and Q^, but 
added the argument y. To get the average momentum trans- 
ferred in a single scattering act, one needs to divide these 
expression by the total cross-section ~){x,y, 77). The cj com- 
ponent of the transferred momentum for isotropic electrons 
is shown in Figure [5^. As shown by NP94, the low-energy 
(Thomson) limit is given by x(l + 2p^/3). The angular de- 
pendent corrections arising due to the dipole and quadrupole 
term in the electron distribution are shown in Figures ^ and 



15};, respectively. While the component perpendicular to co is 
zero for isotropic electrons, a substantial momentum compo- 
nent arises in the anisotropic case. For a large linear term of 
the electron distribution with fi I fo — 1, the momentum trans- 
ferred in that direction is shown in Figure [6^. Similar results 
in case of the quadrupole term with /2//0 = 1, are shown in 
Figure|6j). In the Thomson Umit, we get (see Appendix ID]) 



63 soix, y,ri)^x 



1 + |/-T'7y6^(4r' + 1) 
J Jo 



fi 4 , 

Jo Ij 



Qi so(x,y,T]) =X7 



fo3\ 5j 5/0'^ 



(95) 



(96) 



4. REDISTRIBUTION FUNCTIONS FOR ANISOTROPIC ELECTRONS 

We would like to reduce the expression for the redistribu- 
tion function (|6]l to a form suitable for calculations. For the 
electron distribution of the form this function should de- 
pend on the energies of incoming and scattered photons xi 
and X, the corresponding (cosines of) polar angles 771 and 77 
as well as the difference in azimuth (p - (p\ (or cosine of the 
scattering angle fj). 

4.1. Integration over electron directions 

The three-dimensional integral over p in Equation ^ dis- 
appears due to the 5-function. For further simplifications we 
can also use the identity 



5(y\ + xi - 7 - x) = y5(x ■ (£^ -H Xj) - £j ■ . 



(97) 



At this stage, we drop subscript 1 with the electron quantities 
and get 



— r 



d^p 



6(T)Uy,Ti,)F. (98) 



where 



r = y{xi - x) - p(x\(jJ\ - xcj) ■ il - q. 



(99) 



The angular integrals in Equation ( |98] l need the introduc- 
tion of a suitable coordinate system. Often the polar axis is 
taken along the direction of the scattered photon oi (see e.g. 
iNagirner & Poutanen 1993, 1994). However, the easiest and 
the most transparent way, is to choose the polar axis along the 
direction of the transferred moment um as was propose d by 
lAharonian & AtovanI ( 119811) (see also lPrasad etal.ll 19861) 



n = {x\a)\ - xai) IQ, 



(100) 



where 

i2 



Q —{x\a)\—xco) -X H-Xj -2xxi/i = (x-xi) +2q. (101) 
With this definition we get: 



co?,K = n- M — {x\p - x) I Q, sin/c - xi - fi^/Q (102) 
and 



cos a = n ■ I2 — (x\rji — xrj) / Q. 



(103) 



Thus one of the integration variables becomes cos 9 - Sl-n 
and another is azimuth O. The redistribution function 
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) = Jpdrjd COS eJdO /e(r, ?7e) F 5{T\ ( 104) 



then can be written as 

oo 1 

R{x\ — > x)-- 

1 -1 
where now 

V ^y{x]-x)-q- pQcosO. (105) 
Integrating first over cos 9 using the 5-function we get 



R{x\ — > x) 



3 1 

1^1 '^e Jo ' 



da)/e(r,'7e)^', (106) 



where we need to substitute 



cos Q - 



yjxi -x)-q 

pQ 



(107) 



to the expressions for tj^ and F (see below). This yields 

b 



sin 9 - 



(108) 



where 



V~rJp^Q^-[7(xi-x)-q]\ r=|^. (109) 
> 1 -fi 

The lower limit for the integral over y comes from the require- 
ment that I cos 9\ < 1: 

y>y,{x,xi,ii) = +eVl +2/<7). (HO) 

4.2. Integration over the azimuth 

In order to calculate the azimuthal integral in Equa- 
tion d 1061 1 we have to express ^ and ^\ (that enter the ex- 
pression for F) and rf^ in terms of the integration variable O. 
We measure the azimuth 4> from the projection of co onto the 
plane normal to n, so that in this system 

CO - (sin /c, 0, cos k) (HI) 

and the unit vector along the electron momentum is 

Q - (sin cos O, sin 6* sin O, cos 0). (112) 

Thus we can express the angle between the electron momen- 
tum and /3 (see Fig.[TJ through 

= Q, ■ li - COS0 cosa -I- sin0 sin cos(;^' - O), (113) 

where is the azimuth of the vector 1^ in the n frame. We can 
also write 

T] = a> ■ It, - cos K cos ff + sin sin a co^x^ (114) 

and use t his ex pression to obtain cos;^^. Substituting it to 
Equation (II 13b we thus express the electron polar angle rj^ 
in Equation (O through the integration variable 

The kernel F depends on the four-products ^ and ^\ , which 
can be rewritten as 

^i^x{j-pO, ^ = q + ^i, (115) 

where 

^ = il ■ cj - cos9 COSK + sin sin /<• COS O. (116) 
Equation dl 15t then can be transformed to 



q 



^1 = ^(fl'--^'COS(D), ^^-^{d^-bco^^), (117) 



by anisotropic electrons 1 1 

where we defined 

d- - y(x + xi) - x(x - xifi) , 

d+ - y{x + x\} + x\{xi - xfi) - d- + Q^, (118) 

which have the following property: 

{d^_-b^)/Q^^iy-xf + r = at 
{dl-b^)/Q^^(y + xif + r = al. (119) 



Function F in the azimuthal integral in Equation ( 1106b is an 
even function of Therefore the terms in /e containing sin <1> 
give zero contribution. Neglecting these terms we can express 
the azimuthal integral as 



77eFdO= ^fdO, 
Jo Jo 

7jlFd(i)^ I rilFd<:>, 
Jo 



(120) 
(121) 



where 



?7e =cos 9 COS Q- + sin0 sino- cosx cos (£>, (122) 

rjg — cos^ 9 cos^ a + sin^ 9 sin^ a sin^x 
-H 2 sin sin Of cos^ cos a cosx cosO 
-i-sin^0 sin^a cos2;t' cos^O. (123) 

Thus the expansion ^ (with T/e and rjl substituted by rj^ and 
rjl, respectively) is a quadratic function of cos O. Expressing 



cos O - 



2b q 



2b 



(124) 



cos O = — T- — 5 + — — ; — 

b^ (f- 2b^ q 2b^ 

and using the identity ^ = -H ^, we obtain an expansion of 
/e that is symmetric in ^ and ^\ : 



/e(r) = Co + C2 + Cn^. 

q q 



(125) 



The coefficients cq, cy. and cn can be represented in the form: 

(126) 



Co =/o + coi/i + C02/2, 

CS =Cll/l + C12/2, 
Cn = C22/2, 



where the coefficients in front of /o,i,2 can be easily derived 
after lengthy but straightforward calculation: 



coi - 



Cll 



2p — e + fi 
2p(\+li) ' 



2p(l+//)' 

3 1 

" 4p2( 1+^)2 - Z')' + + Z')' + + - 5' 

3 

"^12 = -T^Ti + fi)'^2p - e -I- ei) -I- /l(<i- -h <i+)] , 

4/?^(l -H /z)-^ 
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Here we defined 

e = x(r]i-j]ij.), 61^x1(77-771//), p^yiT] + r]i), 
A = fi^ + r]^ + r]\-2iJj]7]i-l. (128) 
The redistribution function (1106b is then expressed as 

R(xuMi ^ X, «) = ^ J dj [coRo + c^R^ + cnRn] , (129) 

y,(x,xi,ii) 

where we have introduced three functions 



1 r" 

Jo 
1 r 

rrQq Jo 

1 r 

Jo 



^o( 

/?s(x,xi,/i,y) 
7;n(Jc,-^i,A',r) 



d(D, (130) 
(^ + ^i)FdO, (131) 
f^iFdiB. (132) 



Alternatively, we can represent the redistribution function as a 
sum of three terms arising from the corresponding three terms 
in the electron distribution: 

00 

Rixuioi ^ X, 6>) = ^ jdy [foRo + fiRi + fiRi] , (133) 

y,(A-,A-i,/i) 

where 

/?i(x,77;xi,77i;//;7) = coi/?o +c\\Ry. , (134) 
7?2(x, 77; xi , 771 ; y) = C02R0 +cxiRy. + cziRn ■ 
Using the Klein-Nishina cross-section (|2]i in the form 
q^-2q-2 I 1 1\ 1 1 



F = 2 + 



q 



^1 d'^e^e,' 



(135) 



(and remembering that ^ - + q), we see that the integrals 
( |130l l-( fT32] i involve integrals of types 



r*7T r*7r 

Jo Jo 



(136) 



where s = -2, ..,2. The integrals over non-negative powers 
of ^ a nd g i are trivia l. For the n egative powers, u sing Equa- 
tions dl 17b and (II 19b we get (see lNagimer & Poutanen., 19931 
for details): 

■''dO nQ 1 pdd) nQd 



r dO _ TTg 1 

Jo q a- 



jo q a- Jo ai 

and similar equations for ^ which we get by substituting ^, 0+ 
and £/+ for^i, a_ and £/_, respectively. 

After some straightforward algebra we get the expressions 
for R[), and Rn. 

' •''-''-'(^-^Al'-.^'-il (138) 



Ro^- + 



Q (f- \a- a+j q^\al a 

which was obtained bv | Aharonian & AtovanI (Il981h (see also 
iNagirner & Poutanenlll993h . 

and 



Equation (11291 ) (or alternatively equations 113311 and 113411 ) 
together with our comput ed re d istrib ution functions ( 11381 )- 
( 1140b and the coefficients ( I126l )-( ll28b give the full analytical 
solution for the redistribution function describing scattering 
of arbitrary photons from the electron gas which anisotropy 
can be described by Equation (0. 

4.3. Alternative redistribution functions 

An alternative expression for the redistribution function 
R{xi — > x) can be obtained if we compute the moments 

1 r 

R^(x,xi,yL,y) - — cosOFdO 
^2 Jo 



q -2q - 2 I Id 



b\a 



dA b 
a 



2 1 3" + — 

q \fl_ ai 



(141) 



R^^{x,x■^,^Ji,y)^ — cos^OFdd) 
Jo 

_ 1 _^ 2^ / 2\ q^-2q-2 1 (dl 



Q 



q 



q 



b^ 



a_ a+ 



q^b^ \a_ a+i q- \a^_ a; 
The expressions for Ri and R2 then take the form: 
Ri =doiRo + d\\R^ , 

^2 = 1^02^0 + 1^12^0 + d22R(P(p , 

where 

do\ =cos0 cos Of, 



(142) 



(143) 



dii =sin6' sinci; cos;^^ 



(e + ei), 



do2 = - (cos^ 6 cos^ a + sin^ 9 sin^ a sin^;^') - - 



3 

- ^ . 2 

= ^ (2c/^, - t/^i + sin^ e - cos^ «) - (144) 

di2 = 3 cos 6 cos a sin 6 sin a cosx - 3 doi dn, 
3 

d22 = - sin^ 9 sin^ a cos 2^ 

= ^ (2^/11 - t/oi - sin^ 8 + cos^ a) . 

4.4. Approximate redistribution functions 

Approximate forms of Equations (I138b - (ll40b can be ob- 
tained by making certain simplifying assumptions about the 
scattering. For example, in the Thomson regime in the elec- 
tron rest frame the Klein-Nishina kernel F is just 1 + jJ^y As- 
suming further isotropic scattering in that f rame a nd sub stitut- 
ing F by 4/3, we now get for the integrals ( ll30l )-( fT32] i: 



Rqk — , 
32' 

4 

^2 ~ id- + d+) , 

4 / b^ 

Rn ~ — r d-d+ H 

3Q^\ 2 



(145) 
(146) 
(147) 



The expressi on for Rq was derived by 
lArutvunvan & NikogosvanI (0*9 80.) . For the alternative 
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functions ( 1141b . ( 11421 ). we then have 



These then give 

Ri » doiRo - cos cos a Rq, 
1 



(148) 



(149) 



Ri ~{do2+ 2^^22 1 ^0 = P2(cos 6) P2(cos a) Ro, (150) 



with cos and cos a given by Equations ( 1107b and ( 1103b . re- 
spectively. The approximate expressions are better than 50 
per cent accurate in the Thomson regime for jciy < 0.1 at all 
scattered photon energies. 

4.5. Relation to the mean powers of photon energies 

The relation between the redistribution function averaged 
over any electron distribution and the mean powers of photon 
energies follows directly from their definitions (|6]l, dill and 
(O: 



x\ so(x, rf) 



d^wi R{x xi). (151) 



This relation is vali d for any electron distribution. Compar- 
ing Equations ( 1133b and ( |49] ). we get a relation between the 
functions depending on the electron energy: 



Itlnypx-i^^ 

I x^^^dx\\Rk{xi,riux,ri;n;y)d?'cjJi, (152) 

Jx-(x,y) J 



where rji - t]ij+ -y/T^^^ cos <1) and Rq depends only 

on the scattering angle /u, but not rj,?]]. The integrals over 
the solid angle can be represented as the integrals over dfi 
and dO, where <i> € [0, 2;7r] and the limits on fi, yUm(xi,x, y) 
and fi^(xi,x,y), are given by Equations ( IF6l )-( lF8b with the 
arguments x and xi reversed. Using Equations dSTl l and (|84] |. 
we also get 



A\,(x,y)Pkiv) ^ 



Af,ix,y)Pl{ji) 



32nypx^ 
3 

32nypx- 



JA'^xiJ (l-^)Rk(fLOu (153) 

4' 



'xjdxi IJl - fi^cos<iRkd^a>i. 



(154) 

In order to check the accuracy of o ur der i vation s we com- 
pared the left hand sides of Equations ( I152b -( ll54b to the right 
hand sides, where the integrals were performed numerically 
and obtained consistent results. 

5. APPLICATIONS 

5.1. Examples of redistribution functions 

Now we demonstrate the properties of the derived redis- 
tribution functions. We consider a volume filled by electrons 
with the angular distribution given by Equation (O. The emis- 
sivity in a direction a> at energy x can be obtained from the ra- 
diative transfer equation (|5]l and is given by the integral over 
the redistribution function 



e(x) - crjNe X 



Jo x^ J 



d LOi I{xi) R(xi,a}i X, (jj), 

(155) 
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Fig. 7. — Redistribution functions for anisotropic electrons at a given scat- 
tering angle. The incident photon energy is ,ri = 10"^ and the electron 
momentum p = 0.1. The upper panel shows the photon (number) emis- 
sivity for isotropic electrons. The solid, dotted, dashed, dot-dashed and 
dot-dot-dashed curves correspond to the cosine of scattering angle /j = 
-2/3,-1/3,0, 1/3,2/3, respectively. The lower left panels show the ratio 
R[/Ro, while the right panels show R2/R0 as a function of the ratio of the 
scattered to the incident photon energies. The three row of panels con'e- 
sponds to the different observer directions = 1, 0.5, 0. 

where lixi) - 2mJ^msC^lhf'x\n{xi) is the specific intensity 
of the incident radiation normalized to the photon density as 

A^ph = ^ r r — ■ (156) 

meC^ J J X 

Let us consider mono-energetic (with energy y) electron 
distribution (fTOl l. Consider also a monochromatic source of 
isotropic seed photons at energy xi with total photon num- 
ber density A^ph. According to Equation ( 11331 ) we can write 
the emissivity at an observer direction 7/ for a given scattering 
angle as 



3 X I 

e(x,r],fi) - -r—m^c' o-jNeN^h 

3271 xi py 



f\- fi- 
Ro + -^Ri + -fRi 

JO Jo 



(157) 
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Fig. 8. — Same as Figurel?] but for p = 1. Note, that here the axes are in 
logarithmic units. 

which is related to the scattering angle-averaged emissivity as 
€(x) - J J €(x, T], ;u)dyu, and where (for k - 1,2) 

1 r^" 

Rk{x,ri;x\;n;y) ^ — \ d<S> Rk(x,T];xi,i]uiu-y) (158) 
27r Jo 

and T]i = Tjn + -y/r^^^ -y/r^^/? cos <1>. These functions obvi- 
ously possess symmetry properties: 

Ri{x, -Tj; xi -Ri(x, tj; xi ; yu; y), (159) 

Riix, -i];xi; n; y) = Riix, tj; xi; (160) 

We compute separately the emissivities resulting from three 
terms in t he ele ctron distribution, i.e. functions Ro,Ru Ri (see 
Equation 11571 ). and show in Figures l7l-[T0l the function Rq 
multiplied by x (i.e. quantity proportional to the photon num- 
ber emissivity) for better visibility as well as the ratios Ri/Ro 
and RiIRi). The main behavior of t he fu nctions can be easily 
understood using formulae ( ll49b -( fT50l l derived in Thomson 
limit and isotropic scattering approximation. Averaging them 

over the azimuth and using relation Pkiili) - Pk{ii)Pk(lA, we 
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Tj COS 6 , 



R2 XjP20u) - 2xxijU + jc 
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Same as Figure[8] but for p — 10. 



(161) 



PiiT]) PiicosO) . (162) 



These approximate expressions become extremely accurate 
for high p (i.e. accuracy is about 10"-' at p = 100). 

For a small electron momentum p - 0.1 and low photon 
energies xi = 10"^, the exact redistribution functions are 
shown in Figure |7] In this regime, scattering is nearly co- 
herent with the scattered photon energies bounded by (see 
Equation (IF2l i) x^/xi a: 1 ± p yj2{\ - p). In this regime, 
<s <K X, xi and cos 9 ~ (xi-x)/ pQ is a nearly linear 
function of x/xi, because g/xi ^ ^J2(l - p). For ju not too 
close to 1, the azimuth averaging of cos a gives -// -^(1 - p)l2 

and thus R\/R[) ~ (1 - x/xi)Ti/(2p). For 77 = 0, the func- 
tion is always zero, because of the symmetry. Similarly, the 
nearly quadratic dependence of R2/R0 on energy results from 
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Fig. 10. — Same as Figure[8] but for p = 100. 



x/x-f 



the COS- term, while at » 1/3 the function becomes more 
complicated because of the cancellation in the P2(cos a) term. 

In the opposite limit of the relati vistic electrons (see Fig- 
ures |9] and [TOt . the approximation ( I145l l for the function 
works fine up to xiy < 0.1, while as said above the ratios 
R\/R q and RiJRo are very close to those given by Equations 
(II6II 1 and ( 11621 ) for any photon and electron energies. At 
small scattered photon energies x <^ xi, xR() oc xjxi, and 

^1 IRi) ^ rill cos 6, R2/R0 ~ PiiiDPiiiijPiicos. 0), (163) 

with cos K 1 - xjxi. At high scattered photon energies x » 
jci, the photons are scattered at large angles in the electron 
rest frame and therefore they are beamed in the direction of 
the incoming electrons. In that case, the angular distribution 
of the scattered photons resemble that of the electrons. In 
this regime xRq oc const, and R\IRq ~ ri and R2/R0 ~ PiiTf), 
which gives the flat dependences clearly seen in Figures |9]and 
[TOl and e(_x) oc /^(y, t])x/xi. 



5.2. Sunyaev-Zeldovich effect 

Let consider a cloud of isotropic Maxwellian electrons 
of temperature = kT^lm^c^, which moves with ve- 
locity cySb (corresponding Lorentz factor Fb) through the 
isotropic cosmic microwave background of temperature 
©cmb = kTcmb/meC^. We comp ute the thermal and 
kinematic Sunyaev-Zeldovich eff'e cts jZeldovich & SunvaevI 
Il969f: 'Sunv aev & Zeldovichlll972h . i.e. the spectrum of the 
scattered radiation (and resulting deviations from the black 
body) as a function of and the angle between the line of 
sight and the direction of motion. 

One approach would be to make a Lorentz transformation 
of the incident radiation to the comoving frame, compute the 
Compton scattered radiation using the kernel corresponding 
to isotropic electron distribution, and then to Lorentz trans- 
form it back to the observer frame. Another way is to compute 
the electron distribution in the observer frame, approximate it 
by the expansion ([8]) and compute directly the Compton scat- 
tered radiation in the observer frame. The second approach 
might be favorable from numerical point of view if the ob- 
ject velocity is variable in space and/or time, as allows to pre- 
compute redistribution functions at a fixed grid of angles and 
photon energies. 

5.2.1. Scattering in the comoving frame 

Let us first compute the scattered radiation by a standard 
method considering scattering in the comoving frame. The 
relativistic Maxwellian distribution of electrons in the comov- 
ing frame (quantities with primes) is given by 



/e'(/) = K 



, exp(-y70) 
47r0/:2(l/0)' 



(164) 



where K2 is the modified Bessel function and N'^ is the elec- 
tron density in that frame. The incident black body radiation 
occupation number is 



1 



exp(jf,) - 1 ' 



(165) 



where jc, - jc/0cmb = hv/l<:Tcmh- From the radiative transfer 
equation in the limit of small optical depth, we get the 
correction to the black body spectrum: 

An(x, 77) - n(x,ri)-n},\,(x) - -T'Y'so(x')n^h(x)+S(x,ri), (166) 

where tj is the Lorentz invariant optical depth for Thomson 
scattering, 

S{x, rf) - Tj— J x\dx[ J (fcij[ R"^°{x\ — > x')n\,b(xi) 

(167) 

is the source function, and we used here the fact that the pho- 
ton occupation number is Lorentz invariant. The energy trans- 
formation is given by Doppler shift x - x"D and x\ - x\'D\ 
with the Doppler factors 



1 



rb(i-j8b'7) 



Di ^rbd+ySb^?;). 



(168) 



The relation between the angles is given by the aberration for- 
mula: 

i = f^. (169) 

We note here that io is equal to unity with high accuracy, be- 
cause scattering is in deep Thomson regime. The calculation 
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Fig. 11. — Deviation from the black body spectrum of cosmic microwave background radiation with T^mb = 2.7 K resulting from the Compton scattering in 
a moving cloud of isotropic hot electrons (thermal and kine matic Suny aev-Zeldovich elfects). The electron temperature is = 0.03 and the cloud velocity 
ySb = 0.01. The solid curves are computed using Equations )166) - fT67l . considering scattering in the cloud fr ame, where electrons are isotropic. The dashed 
curves show the results using the formalism developed in this paper for anisotropic electrons, given by Eq uation s {TtD-GH). The dotted curves correspond to the 
semi-analytical approximation of the angle-averaged redistribution function given by Equations UTSt - fTTS^ . The three different methods give nearly identical 
results. The thre e curves from bottom to the top correspond to the three viewing angles with = -0.98, 0, +0.98. The dash-dotted curves show analytical 
approximation of lSazonov & SunvaevI (1998), which includes terms up to second order in /Sb ™d 0, as well as a cross-term It works reasonably well up to 
the temperatures < 0.02, but fails at higher temperatures in Wien tail. 



of the redistribution function involves numerical integra- 
tion over the Maxwellian distribution (see Equation 113311 ; 
note that /, = /z = 0) /o = f^(y) = f^(p)/K given by 
Equation ( 11641 ). Thus the source function ( 11671 ) involves 4- 
dimensional integral to be taken numerically, which is rather 
time-consuming. 



5.2.2. Scattering in the external frame 

We can also compute the same effect directly in the external 
frame. The electron Lorentz factor in the comoving frame is 
related to the electron four-momentum in the external frame 
as 



where y6, = Phi® and we expanded the expression up to the 
second order in ySb- The electron density in that frame is: 

A^e = J /e(p)dV = TbA^e- (172) 

The corresponding terms /j. of the electron distribution can 
be obtained from Equation ( I171l i noting that /e(y, T/e) - 
feipM^hN'^. The change to the occupation number is: 



An(x, rf) - -T-YSoix, T])nbh(x) + S (x, rf). 



(173) 



The scattering cross-section is given by Equation (|26] | and in 
Thomson limit is just 's(){x, rj) ^ \ - fSbrj. The source function 
is now 



r' = rb(r - pPhrie), 



(170) 



S {x, rf) - tt 



where //e is the cosine of the angle between the electron mo- 
mentum and the direction of cloud motion. Because the distri- 
bution function is Lorentz invariant, we easily get the electron 
distribution in the external frame: 

Up) = = ^e^^^^j^ exp(/, Fb/^b l^J&) (171) 

exp(-2V0) [. yS?, 2 . n^^« ^P^P\ t ^ 



X Jo 



xi rn,h(xi)dxi I d (jji R(xi ^ x), (174) 



where the redistribution function R given by Equation ( 1133b is 
averaged over directions of incident photons, but still depends 
on the scattered photon direction 77. This form of the s ource 
function is more favorable compared to Equation (11671 ) from 
numerical point of view, as it can be tabulated in advance at a 
given grid of photon energies and angles. Computed directly 
it still involves numerical calculations of 4-dimensional inte- 
grals. 
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5.2.3. Isotropic scattering in Thomson regime in the electron rest 

frame 

In Thomson limit (as in the case of Sunyaev-Zeldovich ef- 
fect), the calculations in the external frame can be dramati- 
cally simplified. We ca n use t he azi muth ally averaged approx- 
imate expression ( 1145b . ( 1161b . and ( 11621 ) for the redistribution 
functions: 



d R{x\ x)-2n 



1 ^ 



Ro + fiRi +/2-R2]: 



y.(x,xi,tt) 

(175) 

Interestingly, Rq does not depend on y and in expressions for 
Ri and R2 it comes only through cos 6 x (xi - x)j/Qp (be- 
cause q <s x,xi, s ee eg . II107I ). For the electron distribution 
given by Equation ( |171b , the integrals over y thus can be taken 
analytically: 



dyfoRo^C^ 

i b 



/ 

J dy/i/?i=C/?b/7 

/ 



XlyU — X X\ — X 



Q 



(176) 



Pi 



dyfiRi = C — P2(n) g2 



where the proportionality coefficient C - Rq exp(-y»/0)/ 
\4nK2(l /&)]■ T he zerot h order term in wa s deriv ed by 
iPoutanenI (11994 . see also lPoutanen & Sve ns sonI (119961) . 

Evaluation of the source function (1174b now involves only 
two numerical integrations over the photon energy xi and co- 
sine of the scattering angle /i, reducing the computational time 
by 2-3 orders of magnitude. 

For all three methods we numerically compute the correc- 
tion function for the black body intensity 



1 ^ 

A/(x) - — X, An(x, 77), 



tt 



(177) 



and compare the results of calculations in Figure[TT] The three 
different methods give nearly identical results. 

6. CONCLUSIONS 

We have developed the exact analytical theory of Comp- 
ton scattering by anisotropic distribution of electrons that can 
be represented by a second order polynomial over cosine of 
some angle (dipole and quadrupole anisotropy). For the total 



cross-section, we reduce the 9-dimensional integral to a sin- 
gle integral over the electron energy. Analogous expressions 
have been derived for the mean energy of the scattered pho- 
tons and its dispersion. We also obtained analytical expres- 
sions for the radiation pressure force acting on the electron 
gas. These moments can be used for analytical estimations as 
well as for the numerical solutions of the kinetic equations in 
the Fo kker-Planck approximation (see e.g. lVurm & PoutanenI 
I2OO9I) . 

Furthermore, the expression for the redistribution func- 
tion describing angle-dependent Compton scattering by 
anisotropic electrons is reduced to a single integral over the 
electron energy. Exact analytical formulae valid for any pho- 
ton and electron energy are derived in the case of monoen- 
ergetic electrons. We have also derived approximate expres- 
sions for the redistribution function, assuming isotropic scat- 
tering in the electron rest frame, which are very accurate in 
the case of relativistic electrons interacting with soft photons 
in Thomson regime. 

We applied the developed formalism to the accurate cal- 
culations of the thermal and kinematic Sunyaev-Zeldovich 
effects for arbitrary electron distributions. A very similar 
problem arises in outflowing coronae around accreting black 
holes and neutron stars, where the bulk motion causes electron 
anisotropy. Another application could be a computation of the 
radiative transport in the synchrotron self-Compton sources 
with ordered magnetic field, where the electron distribution 
can have strong deviations from the isotropy because of pitch 
angle-dependent cooling. These problems will be considered 
in future pubhcations. 

This work was supported by the CIMO grant TM-06-4630 
and the Academy of Finland grants 122055 and 127512. 
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APPENDIX 
A. FUNCTIONS Sj AND Sj 

All function sj and S j can be expanded to the series, which converge in the region ^ < 1/2. It is easy to show that 

SjiO = 2 aj„ (-2 a\ S j(^) = 2 Aj„ (-2 0", 

where 



n=0 



n=0 
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« + 2 + 
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_ 1 



«(« -I- 5) -I- 



24 
« + 3 



fl2« = ^ («^ -H 9 H- 22 n -H 32) . 



(Al) 



(A2) 
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Using Equations ( l48T l. one can obtain the expressions for the coefficients A through a: 



A\n-2{a\„+\ - flOn+l), - 2 {A\„+\ - fll„+l), A^n - 2 (fl2n+l - Oln+l), (A3) 

Ann-2iAi„+\ - A\„+\), Asn - 3 A4,, - 4A3„, A-i„ - A3„ - A4„/2, A6„ = fl2« - 3 A7„. 



The coefficients for « = 0, 1, 2 are presented in Table[T] 

B. AUXILIARY FUNCTIONS ipij AND T,; 

The total cross-section and mean powers of energy of scattered photons are expressed through the functions of one variable 



i+l ■ 1 dx i+\ 



Calculations of functions i/'/y involve integrals of the following types: 



/d..'«ln(1.2.), /d.^, = f 



djc 

ln(l +2x) — , 



(B2) 



where m = -1, 0, 1, 2, 3 and n,l - 1, 2, 3,4. All integrals are elementary except g{^), which is described in details in AppendixICl 
The explicit expressions for the functions i/^y are the following: 
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= 167? - 4f ' + 5 ^ - 3 /f + ^ ( 1 + 4^ + 2^2) ^3 j ^ 
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32 



11 



18^ - -^f + lOf- -12^+— 1^+^(1+7^+ Uf) 4 



where = ln(l + 2^) and 7?^ = 1/(1 + 2^). 
The exphcit expressions for can be obtained using definitions ( l48b for 5 j: 



*P11 


= 2 (foo 


-<Aoi)/^, 


^13 = 


2(^01 - i/'02)/^, 




2(2iA-ii 


- tff-m - ^-n)l^. 




= l/'l2 - 


3'Pi3 + 3'Fi4/2, 


^21 


= 3(<Aio-«Aii)/2^, 


^22 = 


3(<Aii- 


'Pii)/2f, 


'Pas 


= 3(iAii 


- <Al2)/2^, 


*I'24 


= 3(^ii-^i3)/2^, 


*I'25 = 


3*I'24- 


4'P23, 




= lA22 - 


31'23 + 3>F24/2, 


*P31 


= 4(1^20 -<A2i)/3^, 


*I'32 = 


4(tA2i - 


'P2l)/3f, 




= 4(1/^21 


- -A22)/3^, 


*I'34 


= 4(»I'21-*i'23)/3^, 


"^35 = 


3^34- 


4'P33, 


>P37 


= ^33- 


*P34/2, 


*P36 = 


i/'32 - 3*^37, 


^>4l = 


5 (lp30 - 


-'A3i)/4^, 


*P42 


= 5((A3i 


-*P3i)/4^, 


*P43 


= 5 ((^31 - <A32)/4^, 


*P44 = 


5 CP31 - 


*P33)/4^, 


*P45 


= 3>P44 


-4^43, 


'P47 -- 


= >P43 - *P44/2, 


*P51 = 


6 (ij/40 - 


<A4l)/5^, 


*P54 


= 6(>P4i 


-*P43)/5^, 


*P57 


= 6(lA4I - 1/'42)/5^- 


T54/2. 







(B4) 

In these formulae the argument^ is omitted. For complete evaluation of these functions we need to compute 18 diff'erent functions 
ijjij given above. 

To preve nt th e loss of accuracy if ^ i s ver y small, we can use the series expansions (see NP94) that directly follow from the 
definitions (IBU and Taylor expansions (lAlk 



.A,7(^)^ya,„(-2^)" . ' i>-i, .A-u(^) = y«,„+i(-2^r— 

i + n+ I ^ n + \ 

«=0 n=0 



/7=0 



(B5) 
(B6) 



with Ujn and Ay„ given by Equations (IA2l i and ( IA3l l. respectively. 



C. AUXILIARY FUNCTION 



Calculations of function xj/ij from Appendix IB] involve integral 



g(0 



f 

Jo 



ln(l +2x) — . 



(CI) 



We repeat here for completeness the method of calculations of this integral from NP94. It is possible to write a relation between 
the values of this function on ^ < 1/2 and ^ > 1/2. Let us define for that the auxiliary function for ^ < 1 



It can be presented by series 



g,(0 = 8(^/2) 

\n-\ ? 



ln(l +x) 



dx 



(C2) 



Vx-ir'^ if^<^.<i, 

n-l 

^+ln2 1n^+y^l^y^ if^.<^<l. 



(C3) 



j(:=0 



As ^. we can take 0.8 - 0.9. Then 



(8^(20 
1? , 



if < ^ < 1 /2, 
if^= 1/2, 



(C4) 



+ -ln^(2f)-g.(l/2f) if^> 1/2. 
o / 
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D. ASYMPTOTIC EXPANSIONS OF FUNCTIONS Xj„ AND Aj„ IN THOMSON LIMIT 

Using Taylor expansion (IB5b of functions i^on for small arguments, it is easy to get an expansion of functions ;kt)„ in Thomson 
limit xy <K 1 : 



Xo, 



i„(x,7) = r" ^ (-2xy)' aoiK„ 



where 



2p(l + 1) 



k=0 



i2k + iy.(i-2ky. 



_ (l+ygy y 1 

l+l Zj [yd + 



i:=0 



[yd 



and int(jic) is the integer part of x. A few first functions are 



^1 = 1, ^2 = 1 + jP\ K3 = 1 +P\ K4^ 1+2/3^ + ^5 = 1 + 



(Dl) 



(D2) 



(D3) 



The first three terms of the expansion (ID lb are as follows: 



XOn(x,y) ~ / 



26 



K[+n - 2xyK2+„ + — (xy) K3+„ 



(D4) 



Function Aqo coincides with;^foo^ and functions Aoi and A02 can be obtained using definitions (IZTT l and expansion (IDlb . For Aqi, 
we get 



Aoi =-/3^(-2x7)'ao/^; 



where 



/=() 



mt[(/-l)/2] , ^ 



P 



(2^ + 3)!'^ (l-l-2k)\' 



78, ,02' 



1 - 4-xy + — (xy) + — 



(D5) 



(D6) 



4 4 



6 = ^2 = -, ^3 = l + 7y8', ^4 = - + 7y6\ ^5 = " + 2/3^ + -yS^, ^6 = 2 + 4^^ ^ -yS^. 



3 ' 3 
Respectively for A02, we have 



3 5 



A02 J(-2x7)'flo/A/+i * -^/^^(^r) - y -^rj , 



where 



A; = - 



2yS2 



/=0 



ki-2km +K1+2 
3 — Kl 



int(//2) 



/! 



2/t 



(Z-2A:)!(2A:)! (2A: + 1)(2A: + 3) 



2(k-\) 



k=\ 

4 4 



2 2 T T 

Ai=0, A2 = — , A3 = -, A4 = - + — fll 
^ 15 5 5 35^ 



Similarly, for functions using expansions < IB6b we get: 

CO "CO 

Xin(x,y) = y"^' ^ (-2xy)' (yAi/ + XA2;) /c„+/+2 = r" r^/^2+n + ^ (-2x7)' ^y^ Ai//f„+/+2 y^/fn+z+ij 



(D7) 



(D8) 



/=() 



(D9) 



Functions Ai„ can then be obtained using definitions 

4 2 xy 



Aio =;rio * 1 + -p' - ^ (42y2 - 27 - 2^^) , 



An^-P Z(-2xy)' (yAi, + x A2;) ^ 1 3^' " y [21r'(l +P'/5) - 4]^ , 



2 r 3 

A12 =;'yS^(-2xyy(rAi, + XA2/) A,+2 ^ yjS' r' - ^ (2ir' - 2) 

/=0 13 [ 3 



(DIO) 
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/=0 

[7 



r 1 



(DID 



where we kept only the zeroth term in xy of the series. Expansions for A2„ can then be obtained using definitions ( l5Ql t: 

2 

Aio^Xio ^ 1 + —/J^ (217^+4), 

°° [2 
A21 = -jS ^(-2^7)' [(/A4, - r'Av,) f,+3 - r'As, ^,+2 + Ae, ^z+ij ^-(3 1 + —p" (63r' + 34) 

/=() L 

A22 ^(-2xr)' [(/A4/ - r'A7,) A/+3 - r'^s/ A/+2 + As, A/+i] * (42/ - 1 1) . 



(D12) 



/=() 



Let us n ow discuss the properties of functions The series expansion can be easily obtained from the definition ( 189] ) and 
series ( |B6t : 



i;(x,r) = £(-2xr)'Ai,^„+,+3. 



(D13) 



For A* , we get: 



00 

=-rpX(-2xr)'A„^,+3 ^ -ypii +Ii'l5), 

1=0 



The series expansion for functions j^^j-^^ are: 



xti(X' r) = — ^^(-2-J:y)'Ai, (/f„+/+2 - 2y^ /f„+/+3 + 7^k„+i+4) = ^/^'p ^(-2xr)'Ai, [0^A„+i+2 - k„+i+2) . 



(D14) 



(D15) 



For Af. we get: 



A^i = ^;ri,--^r/'(i+y6V5), 



11°° 2 

E. ELIMINATING CANCELLATIONS IN REDISTRIBUTION FUNCTIONS 



(D16) 



If formulae (I138I I. ( |139l l and ( 1140b are used as they stand, numerical cancellations appear at certain regions of parameter 
space. For example if x and xi are small, the quantities a_ and a+, l/fl_ and l/a+, are close to each other Also a combination 
containing a sum of d-/a^_ and d+la\ minus double the difference l/fl_ and has a cancellation. Therefore it is useful to 
rewrite the expressions in a form not containing those cancellations. The cancellations appearing in ( 1138b were dealt with in 
iNagirner & Poutanenl(ll993h . Defining 



(x + xi)(2y + xi — x) 
u — a+ - 0- — , V — 0-0+, 



fl_ + 0+ 



they got 



2 u 

Q V 



1 \ + u 

q 



Using definitions dETb . we get from ( 11391 ) and ( 11401 ) 

= (o- + a+) 



23 V 



T— 5-5 + "-TT- 

2^^v' q'^y^ 



2m 1/^ 2\ {u'-Q'){u' + 3v) 



(El) 



(E2) 



(E3) 
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= i + 4v) + 2b^] + -L (l - 1) + (E4) 

Another loss of accuracy o ccurs in - term, when y is close to jt{x,x\,ix). We can use the following formulae 
(iNagirner & Poutanenll993[) : 

- -IrqCDu, Du - iy + x\ - x + j»)(y - y^), C - 2l[y(y + xi - x) + r + xxi^ + v]. (E5) 

F. BOUNDARIES 

The redistribution functions Rq, Rj_, Rn and 7? i , are defined within the inte rval o f photon and electron energies and scattering 
angles satisfying the relation | cos 01 < 1, where cos0 is given by Equation ( 1107b . These limits were discussed in NP94, but 
we repeat them here for completeness. For fixed photon energies and scattering angle, we already got the limits on the electron 
energies given by Equation (II lOl i. y > y,(x,xi,fi). If we are interested in the interval of scattered photon energies for the fixed 
xi,y and /i, we have then < x < x^, where 

±. ^ fi + y(y + xi)(l - ^i) + p (1 - ^i) a+ 
X (xi,r,^) = xi — — ^ — 

If the energy of the scattered photon x is fixed, the initial photon X[ lies in the interval 

x7<xi<x| ifO<x(l-ij)<y-p, 
Xi > x^ if y - p < x(l - jj) < y + p. 



xt(x,y,,^)^x , , _ (F3) 



where 

ji + yiy-x)il-n)+p(l-iJ.) a - 
1 -27x(l -fi)+x^ (1 -pf 

In Equations ( IFll i and ( IF3l l. the quantities a± are defined by Equations dl 19l l. If |x - xi| < 2xxi, the quantity yt(x,x\,p) as a 
function of p has a minimum 

rmin = 1 + (X - Xi + |X - Xi |)/2 (F4) 

at yu = /imin = 1 - |x - xil/xxi, while in the opposite case, |x - xi| > 2xxi, the function is monotonic with the minimum reached 
at the boundary p - -I (see Fig. fT2l l. Correspondingly, the limits of variations of p depend on the photon energies x, x\ and the 
electron energy y and are given by 

Pm<H< H+, (F5) 



where 



{-1 if |x - xil > 2 xxi, 
-1 if |x - xil < 2 xxi and y > yt,(x, xi, -1), 
//_ if |x - xil < 2xxi and y < y»(x, xi, -1), 

A^-(x,xi,r)=l-— , (F6) 

XXl 

p+(x, xi, y) = 1 = 1 , 

XXl XXl q+ 

and 

7t(x, xi, -1) = l^x - xi + (x + xi) yjl + 1/xxij /2, (F7) 



q± ^ +y{xi-x) + p ^jiy + xy -x)2 - 1. (F8) 

For the angle -averaged redistribution function, the lower limit on the electron energy is: 

( . ^\y*ix,xi,-l) if |x- xil > 2xxi, 
7* x,xi <y^^^^ if |x-xi| < 2xxi. 

The limits of variation of the scattered photon energy x as a function of incident photon energy x\ and y can be found by inverting 
Equation ( |F9] l. We obtain 

x"(xi,y) < X < Xm(xi,7), (FIO) 
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Fig. 12. — The dependence of the function •y,(;c, xi,/!) on ji. The energy of the incident photon is .vi = 2. The solid curves are for x = 3 and the dashed curves 
are for x = I (both cases correspond to |x - a'i | < 2xxi ), and the dotted curve is for x = 0.3 (where |x - xi | > 2xxi ). 



where 



•y + xi-lifl<y<l + 

2jc2 

x+(jci,r) if 1 + - — ^ 

1 — ZXl 

,y + xi - 1 if xi > 1 /2, 



2x] 
1 - 2x1 



and xi < 1 /2, 



< y and xi < 1 /2, 



x='(xi,7) = xi (y + p)/{y + p +2xi). 



(Fll) 
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